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Abstract: The electron localization in the dissociation of the asymmetric 
charged molecular ion HeH^+ exposed to an intense few-cycle laser pulse 
is studied by solving numerically the 3D time-dependent Schrodinger 
equation. By varying the carrier-envelope phase (CEP) and the intensity 
of the pulse, the upward shift of the localization probability and the 
suppression of the dissociation channel He^+H-H are observed . Our analysis 
shows that the phenomenon is attributed to the asymmetric structure 
of the molecule as well as the recollision-assistant field-induced ioniza- 
tion of the electron wave packets localized on in the trailing of the pulse. 
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With the development of laser technology, the molecular dynamics in strong field is progres- 
sively understood in depth by scientists. Many of the processes in terms of strong-field molec- 
ular interactions have been discovered and interpreted in experimental and theoretical studies 
on H2 and its isotopes, such as bond softening, vibrational trapping, enhanced ionization and 
so on |[l]. 

Recently, the electron localization in dissociating molecules has become an interesting sub- 
ject for molecular research ||2l|3]|4|. The control of electron localization in molecules on the 
attosecond time-scale is crucial for steering many physical and chemical processes fT, "^l. In 
a landmark experiment, Kling et al. [7 J reported the probability to control the direction of the 
D+ ejection (and hence the localization of the electron) in the dissociating by varying the 
CEP of the pulse. This constitutes the first realization of direct light-field control of a chemical 
reaction via the steering of electronic motion, and it has stimulated continuous interest to study 
the electron localization in dissociating small molecules ||g[g[Tl[n][T2l[T3l[T4l[T5l[T6l[T7l[T8l 
that include two nuclei with the identical charge. More recently, Znakovskaya et al. |fT9l have 



reported an experimental demonstration of electron localization in the dissociation of a het- 
eronuclear molecule, CO, and the control of the directional emission of C+ and fragments 
has been achieved in their work. 

Motivated by understanding more detail of the physical mechanism for electron localization 
on asymmetric charged fragments, we set our efforts on theoretically studying the dissociation 
of a model molecular ion, HeH^+, in an intense few-cycle pulse. As the simplest asymmet- 
ric charged molecule, the cation of HeH has been the subject of many works on laser-induced 
phenomena, such as high-order harmonic generation ||20lETI . enhanced ionization |22|, double 
ionization |231 and photodissociation |24|. The molecular system of HeH^+ consists of three 
charged particles and a bound state of the system is existent |25|. Interacting with a strong 
laser field, HeH^+ exhibits three fragmentation pathways, the dissociation channels He+H-H"*" 
and He^+H-H, where the electron localizes on either of the nuclei, and the ionization channel 
He^^+H^+e . Hereafter, p and a are used to denote the proton and alpha particle He^+, 
respectively. In the present work, the CEP effect on different dissociation channels of HeH^+ 
is studied by solving numerically the 3D time-dependent Schrodinger equation for the interac- 
tion. The numerical results show an interesting phenomenon that the probability of the electron 
localized on p shifts upwards with increasing laser intensity and the dissociation channel a+H 
is suppressed for some CEPs of the pulse. We then qualitatively explain the new findings by 
analyzing the electron localization during the dissociation process. 

For numerical simulation we have used the three-dimensional model for HeH^+ interacting 
with a linearly polarized few-cycle pulse. This model consists of one-dimensional motion of 
the nuclei and two-dimensional motion of the electron lf3l fTTIl . Thus non-Born-Oppenheimer 
effects are included during the interaction. In our model, the molecular ion is assumed to re- 
main aligned with the linearly polarized pulse. Rotation of the molecule is not considered since 
the molecule does not have time to rotate significantly during the dissociation (a few tens of 
femtoseconds in our simulation) as well as the molecule dissociates primarily along the laser 
polarization. Within this model, the time-dependent Schrodinger equation can be written as 
(e = nie = fi — I) 

i^^^{R,z,p;t) = [Ho + V{t)]^'{R,z,p;t), (1) 
where Hq — Tn + Te + Vq is the field-free Hamiltonian with 

r---— (2) 



2^e ^(9z2 ^p2 ' p 



(33 + 3:3 + -3:7)' (3) 



Vo = ^ ^ (4) 

VR^+b ^{z-Zp^ + p^ + a ^{z-za)^+p^- + a 

and y (f) is the electric potential including the laser-molecule interaction. Here, R is the inter- 
nuclear distance, Cp and Ca are the electric charges of the nuclei, and Zp ~ ma/{mp +ma)R 
and Za = ^iiip/ {mp + ma)R are the positions of p and a, respectively, /i = (1 /nip + l/ma)^' 
and /ig — [nip + nia)/ {nip + nia + 1) are the reduced masses, with nip and nia the masses of 
p and a. The soft-core parameters a = 0.01 and b = 0.12 have been chosen so that the model 
of HeH^+ yields the 2pa bound state energy and equilibrium distance of -0.5317 and 3.9 a.u. 
II25I . respectively. The interaction with the laser is treated in the dipole approximation and 
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Fig. 1. The 2pa bound state of HeH^+ and the schematic diagram of two different dis- 
sociation channels for HeH^+: He++p and a+H. The up and down dash hues indicate the 
positions of the p and a, respectively. The color scale is logarithmic. 



length gauge 



y{t) = fj'^'hL^f^R + (1 + ^4^)z] X £ (r) (5) 
mn + mn mn + m„ + l 



with E{t) = £'oexp[— 21n(2)(f/T)^] cos((»f + 0), where £0 is the peak electric field amplitude, 
T is the pulse duration, (O is the central frequency, and is the CEP of the pulse. The electric 
field is chosen to be linearly polarized along the z axis. Then, one can solve Eq.(l) with the 
Crank-Nicolson method and evaluate the time-dependent wave function. 

Before the time evolution of the wave function, an initial state is requisite. However, a stable 
wave function cannot be obtained directly by imaging time propagation of Eq.(l) because the 
ground state of HeH^+ is unstable. This can be understood as that the vibrational wave packets 
are not bounded by the Ls^cJ curve. In our calculation, the first excited state 2pa, which is the 
lowest electronic bound state of HeH^+, is chosen to be the initial state of the system. Within 
the Born-Oppenheimer approximation, the initial wave function is treated as ||27]| 



1'o(/?,z,p) = Zo(W(z,p;^), (6) 

where Xo{R) denotes the first eigenstate of the nuclei bounded by the 2pa potential curve and 
^e{z,P',R) is the first excited state of HeH^+ at each fixed nuclear separation R. In detail, the 
electric wave function ^e{z,P',R) is found by solving the Schrodinger equation of the electronic 
motion at each fixed separation. The potential curve Vn{R) of the 2pa state can be obtained 
by V„{R) = {^'eiz,p■,R)\iTe+Vo)\'i'eiz,p■,R))■ Then the first eigenstate of the nuclei ;Ko(^) is 
solved with imaginary time propagation of the following Schrodinger equation ||271 

i^^X{R.t) = [T„+VniR)]x{R,t). (7) 

The wave function density distribution for the 2pa bound state of HeH^+ is illustrated in Fig.[l] 
In the present simulation, the grid ranges in /^-direction from to 25 a.u., for z from -45 to 
45 a.u. and for p from to 15 a.u. with 500, 300, and 50 points in the three directions, respec- 
tively. The wavelength and duration of the laser pulse are 800 nm and 5.3 fs (FWHM), respec- 
tively, while its CEP and intensity are variable. We have adopted the approach from Ref.||3| 
to analyze the wave function but defined a different configuration space (see Fig. [T]l. The two 
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Fig. 2. CEP dependence of the electron localization probability on the a (blue dots) and p 
(red stars) for four intensities: (a) / = 4.0, (b) / = 5.0, (c) / = 6.0, and (d) / = 7.0 in units of 
lO'^^W/cm^. The black thick curve indicates the reverted sine function. The green arrows 
indicate the suppressions of the probabilities. 

channels of dissociation are defined as ^ ^ > 12 and ■\/ (z — Za)^ + <7, Q.^ : R > 12 and 
y/ {z — Zp)^ + < 5. In the region Q.2 the electron localizes on the a contributing to the disso- 
ciation channel of He++/?, whereas in Q.^, it localizes on the p contributing to a+H. The region 
describing the ionization is not contained in Fig. 1 . We estimate the ionization probability by 
storing the absorbed wave packets reaching the boundary in the z direction. The propagation of 
the wave function has been continued until the probability of the region ili is converged. Dur- 
ing the propagation, we have employed absorbing boundaries using cos '/^-masking functions 
but stored the absorbed contributions as dissociation or ionization, respectively. 

Figure |2]reveals the probabilities of dissociation into He+H-p (PHe+) a+H (Ph) as a func- 
tion of the CEP for different pulse intensities. On the one hand, the strong CEP dependence of 
the dissociation channel of He^+p is obvious, indicating that the control of electron localiza- 
tion on the a can be achieved by adjusting the CEP of the pulse. Moreover, the amplitude of the 
localization probability is larger for stronger pulse. The probability for the dissociation limit of 
He^+p, which associates with the excited electronic state, is small and increases slightly with 
increasing intensity. This is a result of the need to absorb more photons. On the other hand, the 
CEP dependence of Pn is extraordinary. An upward shift of Ph can be seen as the intensity in- 
creases, and the dependence of Pn on the CEP is not a "sine-like" curve but it is suppressed for 
some CEPs, as indicated by the arrows in Fig. [2] The curves of the reverted sine functions are 
also shown in the figure. Here, the sine function is given by f{x) ^ ao + ai sm[{x + a2)(p] and 
the parameters aq, ai and 02 are fitted with the values of P^. To quantify the suppression, we 
define the absolute suppressed probability as the maximal deviation of the dissociation proba- 
bility from the reverted sine function. The suppressed probabihties are 0.003, 0.005, 0.012 and 
0.02 for Fig.|2|a)-(d), respectively, which means that the suppression is deeper for higher laser 
intensity. 

In order to understand the processes of the different dissociation channels during the inter- 
action, we illustrate the model for the dissociation of HeH^+. First of all, the electron density 
distributions of the lowest two states of HeH^+ at each fixed nuclear distance are depicted in 
Fig.[3ja) and (b), respectively. One can see that the electron density concentrates at the a in the 
Isa state while it concentrates at the p in the 2p(7 state. Thus the localization of the electron is 
related to the energy level of the molecule. If the HeH^+ dissociates following the Isa potential 
curve, the electron would locaUze on the a contributing to the He^+p channel. Similarly, if the 
HeH^+ dissociates following the 2po potential curve, the electron would localize on the p con- 
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Fig. 3. The electron density distributions of (a) the ground state and (b) the first excited 
state of HeH^+ at each fixed nuclear distance. The p -dimension has been integrated, (c) A 
simple model for the dissociation of HeH^+. The dashed line illustrates the dressed 2pa 
state of HeH^+ in a static field, and the Gaussians illustrate the wave packets. The color 
scales in (a) and (b) are linear. 

tributing to the a+H channel. Then we illustrate the dissociation process of HeH^+ in Fig.jsjc). 
First, the wave packets are bounded in the 2pa state. When the laser pulse is applied, the elec- 
tric field dresses down the 2pG potential curve f22\ creating vibrational wave packets. Then 
the populations are transfeiTed between the 2pa and IsG states by the external field during the 
interaction. Pait of the wave packets are ionized when the maximum light field is around. With 
Fig. [3ja) and (b) in mind, the dissociating wave packets along the Isa and 2pa curves result 
in the electron localization on the a and p (and hence the dissociation channels of He+H-p and 
a+H), respectively. 

Based on the dissociation model of HeH^+, the simulation results in Fig. [2] can be qualita- 
tively explained as follows. On the one hand, the electron motion is sensitive to the CEP of the 
few-cycle laser pulse. When the CEP varies, the population transfer between the Isa and 2pa 
states during the interaction will be modulated, resulting in the CEP-dependent probabilities 
of the two dissociation channels. On the other hand, according to Ref [22|, the 2pa curve is 
dressed down in both cases of £ > and E <0. Therefore, the 2po state of HeH^+ is always 
dressed down as the field oscillates. The 2pa curve would be dressed even lower if the inten- 
sity of the pulse is higher [22|, leading to more dissociative wave packets on the 2pa potential 
curve. With Fig.[3jb) in mind, this means that more electron density concentrates at the p in the 
dissociation process and thus the probability of dissociation channel a+H shifts upwards when 
the pulse intensity increases. In contrast, the electron localization probability on the a does not 
shift upwards because of the asymmetric distribution of the wave packets in the initial state. In 
the initial bound state (see Fig.l), the wave packets are primarily localized on the p. Thus the 
wave packets which localize on the fragment a during the dissociation have to be transferred 
from the p by the external field. However, due to the oscillation of the pulse, most of the wave 
packets localized on the a are diiven back to the p if the electiic field points to the — z direction, 
as shown in Fig.4. If the intensity increases, more wave packets will be transfeiTed from the p 
to the a when the field points to the +z direction. But meanwhile, the wave packets driven 
from the a back to the p will also increase when the field reverses the direction. Until the field 
becomes weak, a small part of the wave packets are left around the a because of the deep po- 
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Fig. 4. Time evolution of the wave function density integrated over the p- and R- 
dimensions. The intensity of the pulse is 7 x lO'^W/cm^, and the CEPs are 0.5;r and 7t 
for (a) and (b), respectively. The solid and dotted lines illustrate the electric force on the 
electron and the profile of the electric field, respectively. The dashed lines indicate the 
different times when the field reaches its peak value. The color scale is logarithmic. 



tential well and subsequently dissociate with the a. Therefore, the localization probability P^^+ 
predominately depends on the potential well of the a but it does not increase significantly with 
the intensity. 

Next, we focus on the mechanism responsible for the suppression of the dissociation proba- 
bility of a+H. To gain further insight into the interaction between HeH^+ and the laser pulse, 
the temporal evolution of the electron density is depicted in Fig. |4] The pulse intensity is 
7 X lO'^^W/cm^, and the CEPs are Q.5n and n for Fig. ^a) and Fig. ^h), respectively. The 
solid curve indicates the electric force F{t) — —E{t) on the electron. As shown in the picture, 
after the laser pulse is applied, the bound electron density is partially released when the field 
reaches the vicinity of its maximum. One part of the released electron wave packets are driven 
away and ionized. The other part of them return to the nuclei as the electric field changes its 
direction 1281 |29l [30l [31 1 and rescatter with the bound wave packets. Then we show the inter- 
atomic barrier potential and the time evolution of the probability density for R-dimension in 
Fig. 5. In the present model of HeH^+, the height of the barrier crosses the energy of the 2pa 
state at Rq = 7.5 a.u., as shown in Fig.|5ja). According to Fig.|5jb), the time for the dissociative 
wave packets to reach Rq is at about 9.5 fs. After f = 9.5 fs, the electric field becomes very weak 
and the wave packets transfer between the nuclei is blocked by the increasing interatomic bar- 
rier. Then the molecule dissociates with the electron localizing on one of the two dissociating 
nuclei. This is consistent with our dissociation model of HeH^+. 

For the CEP of 0.57r, as shown in Fig. ^a), the returning wave packets at about f = 6 fs 
interfere with the bound wave packets on the p and then the density distribution is modulated. 
One can see that there are bumps in the bound wave packets, indicated by the black arrow, and 
we deduced that part of the bound wave packets are excited to higher states. At the same time, 
the electric field just reaches its peak value. Although the instantaneous strength of the field is 
weak, part of the wave packets localized on the p after excitation are still ionized by the field, 
leading to the week ionization burst marked in the figure. As a result, the probability of the 
electron localization on the p is suppressed. In the case that the CEP equals n in Fig.|4|b), the 
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Fig. 5. (a) The interatomic barrier potential, (b) the time evolution of the probability density 
for R-dimension. The laser parameters are the same as those in Fig.4(a). The thick line in 
(a) indicates the energy of 2pa state at Rq = 1.5 a.u.. The dash line in (b) indicates the time 
when the dissociative wave packets reach R^. The color scale is logarithmic. 



returning time of the wave packets is not the exact time at which the electric field reaches its 
peak value. Although the electric field at f = 5.3 fs in (b) is stronger than that at f = 6 fs in (a), 
the ionization burst is weaker at f = 5.3 fs in (b). Therefore, the suppression of the localization 
probability does not happen for all CEPs. 

On the other hand, the deepest suppressions of the localization probabilities appear to be at 
different CEPs for different intensities in Fig. 2. It has been shown that the dynamics of the 
electronic motion in molecules is affected by both the external laser field and the diffraction 
effects |12|. In Fig. 2, the regions for the suppression are mainly around 0.5 7t for different 
intensities, which are determined by the external field. However, due to the diffraction effects, 
the electronic motion will be modulated for different intensities. Therefore, the positions for 
the deepest suppression are changed for different intensities. In addition, for more intense laser 
pulse, the electric field is stronger and the returning wave packets take more kinetic energy, 
and thus more electron density will be ionized from the dissociating p. This comes to the result 
that the suppression of localization probability is deeper for higher intensity in Fig.|2] One may 
expect that such ionization burst would appear at the other side if the electric field is reversed. 
However, the potential well of the a is deep and the energy level of the Isa state is very low. 
The electron firmly localizes on the a during the dissociation. In contrast, the energy level of the 
2pa state is relatively higher and the potential well of the p is shallow. The electron localized 
on the p would be easily ionized if it is excited by the returning wave packets. Therefore, the 
suppression occurred in the dissociation channel of a+H rather than He^+p. 

Finally, we note that the CEP dependence of the electron localization in the asymmetric 
charged molecular system of HeH^+ is different from that in the previous studies on HD+ 
||3]|4|, in which the upward shift and suppression of the dissociation probabilities have not been 
reported. Moreover, the electron localizations on the p and a are different by changing n of 
the CEP. This is because the mechanism responsible for the electron localization in HeH^+ 
dissociation is different from that in HD+. For HD+, the asymmetric localization on either 
the p or d (D+) is induced by the coherent superposition of dissociative wave packets on the 
gerade and the ungerade states, i.e., the IsGg and 2/?(7„ states. However, the localization in 
HeH^+ dissociation is the result of the dissociative wave packets on the lowest two states which 
have asymmetric density distributions. Compared to another previous experimental work on 
the dissociation of CO+ lfT9ll . the asymmetries of both the dissociation channels, He++p and 
C++0, where the electron is localized on the larger nucleus, shows the similar feature of the 



CEP dependence. There is no comparison for the other channel since the intensity dependence 
of it in asymmetric molecule has not been reported before. 

In conclusion, we have studied the electron locahzation in the dissociation of the asymmetric 
molecular ion HeH^+ by solving numerically the 3D time-dependent Schrodinger equation. The 
upward shift of the localization probability and the suppression of the dissociation have been 
observed. The physical mechanism responsible for the suppression is that the electron wave 
packets which have already locaUzed on the dissociating nucleus is excited by the returning 
wave packets and subsequently ionized by the field at the trailing of the pulse. The suppression 
of dissociation is also related to the localization of the electron as well as its corresponding en- 
ergy level. Generally, the phenomenon of the suppression occurs in the dissociating asymmetric 
molecular ions where the electron at higher energy state is localized on the nucleus with lower 
potential well, such as HeH'^+ and LiH^+. Though the asymmetric molecular ion HeH2+ in the 
present work is a model system, the electron localization on the larger nucleus shows a similar 
CEP dependence as that in CO+. 
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